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ABSTRACT 

In this paper we provide the first consideration of magnetosphere-ionosphere coupling 
at Jupiter-like exoplanets with internal plasma sources such as volcanic moons. We 
estimate the radio power emitted by such systems under the condition of near-rigid 
corotation throughout the closed magnetosphere, in order to examine the behaviour 
of the best candidates for detection with next generation radio telescopes. We thus 
estimate for different stellar X-ray-UV (XUV) luminosity cases the orbital distances 
within which the ionospheric Pedersen conductance would be high enough to main- 
tain near-rigid corotation, and we then consider the magnitudes of the large-scale 
magnetosphere-ionosphere currents flowing within the systems, and the resulting ra- 
dio powers, at such distances. We also examine the effects of two key system pa- 
rameters, i.e. the planetary angular velocity and the plasma mass outflow rate from 
sources internal to the magnetosphere. In all XUV luminosity cases studied, a signif- 
icant number of parameter combinations within an order of magnitude of the jovian 
values are capable of producing emissions observable beyond 1 pc, in most cases re- 
quiring exoplanets orbiting at distances between '^l and 50 AU, and for the higher 
XUV luminosity cases these observable distances can reach beyond ~50 pc for massive, 
rapidly rotating planets. The implication of these results is that the best candidates 
for detection of such internally-generated radio emissions are rapidly rotating Jupiter- 
like exoplanets orbiting stars with high XUV luminosity at orbital distances beyond 
'^l AU, and searching for such emissions may offer a new method of detection of more 
distant-orbiting exoplanets. 

Key words: Planetary systems - planets and satellites: aurorae, magnetic fields, 
detection. 



1 INTRODUCTION 

, In recent years hundreds of exoplanets have been dis- 
■ covered, many of which (~21%) have mass greater than 
' or equal to that of Jupiter and orbital semi-major axes 
"of < 0.1 AU (where 1 AU « 1.5 x lO" m), although a 
significant fraction (~61%) of planets with Jupiter's mass 
or greater have been observed with semi-major axes ^ 1 AU 
(see, e.g., the catalogue at exoplanet.eu). The possibility of 
detection of the auroral radio emissions of 'hot Jupiter'-like 
exoplanets close to their par ent star has been considered 



by a number of authors (e.g. Farrcll, Dcsch & Zarka 


19991; 


Farrell et al.ll2004 


: Zarka et al,,.2001..,2007:.Lazio et al...2004 


Griefimeier et al.l 
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I2OIOI : iReiners fc ChristensenI |201G| ). This interest has 
been sparked in part by the imminent commencement 
of high sensitivity radio observations by next generation 
radio telescopes such as the Low Frequency Array (LO- 
FAR), which will have a detection th reshold of I mjy 
(where 1 Jansky = 10"^*^ W m~^ H z~^) l|Farrell et al.ll2004i : 
IGriefimeier. Zaxka. fc Spreeuwll2007h . Such previous consid- 
eration of hot Jupiter-like exoplanets has assumed that the 
auroral radio emission would be caused by a star-planet 
interaction reminiscent of either the solar wind-Earth 
interaction or the lo- Jupiter interaction. The former is 
mediated primarily via reconnection of the planetary and 
interplanetary magnetic fields at the dayside magnetopause 
(Dungcy 1961). This drives plasma flows within the magne- 
tosphere that generate electric currents flowing between the 
magnetosphere and the resistive ionosphere, the upward 
magnetic field-aligned component of which, associated with 
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downward-precipitating electrons, produces auroral and 
radio emissions. The latter interaction is thought to be 
mainly associated with the generation of Alfven waves in 
the vicinity of the moon lo, caused by its motion through 
the rapidly-ro t ating planetary magnetic field and plasm a 
ijGoert j Il980l : iNeubauej Il98d : ICrarv fc Bagcnall Il997l l. 
Consideration of such processes has led to the extrap- 
olation of a 'Radiometric Bode's Law' relating incident 
magnetic power to output radio power for the case of 
hot Jupiter-like exoplanets orbiting extremely close (at 
typically ~10 stellar radii) of their parent stars. It has been 
concluded that such interaction may gener ate emissions 
at or above the LQFAR detection thresh old (jFarrell et al.l 
12004 iGriefimeier. Zarka fc Spreeuwl |2007|) . 

However, despite the importance placed by previous 
authors on stellar wind-planet and lo- Jupiter type interac- 
tions, significant components of Jupiter's radio emi ssions, 
i.e. t he b-KOM, HOM and non-Io-DAM emissions (jZarkal 
Il998l ). are thought to be generated by the large-scale 
magnetosphere-ionosphere (M-I) coupling current system 
associated with the breakdown of corotation of iogenic 
plasma in Jupiter's middle magnetosph e re, illustrated by 
Fig. [ J (.Hill ,1979l. I2OOII: iPontiud [l997l: ICowlev fc Buncel 
I2OOII : iNichols fc CowlevI 12001 |2004 |2005| ). This process 
generates intense field-aligned electron beams which drive 
the brightest and most significant of Jupiter's UV aurora l 
emission, i.e. the main auroral oval |Grodent et al. I l2003al : 
IClarke et all |2004| : iNichols et all l2009bl ). and, coupled 
with particle mirroring and the absorption of particles 
in the loss cone, excite the cyclotron maser instability in 
the high-latitude low-/? plasma, which gives rise to the 
above radio emissions. Observationally, the UV aurora 
and radio emissions of Jupiter and Saturn have been 
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one another (e.g. 'Gurnett et al. 




2002 


■ 


Prvor et al. 


2005 


: Kurth et all 2005: Clarke et al. 


200s 


Lamv et al. 


200c 


: INichols et all I2OI0I: INichols. Cowlev fc Lamv l2010h. 



lo orbits deep within Jupiter's magnetosphere at ~5.9 Rj 
(where Rj represents the equatorial radius of Jupiter equal 
to 71,373 km), and its volcanoes liberate sulphur and 
oxygen atoms into a torus surroun ding the moon's orbit 



at the rate of •^1000 kg s (e.g. iHill. Dessler fc Goertz 



19831 : IVasvliana3ll983l: iKhurana fc Kivelsonlll993l : iBagenal 



I997I : iDols. Delamere fc Bagenall I2OO8I '). These atoms are 



ionised by electron impact ionisation and thus become 
sensitive to the rotating planetary magnetic field, such 
that the newly-created plasma is picked up to corotation 
velocity. The picked-up plasma is centrifugally unstable 
and diffuses radially away from the planet, probably via 
flux-tube interchang e motions llSiscoe fc Summers! Il98ll: 
Pontius fc ffiii 1 19821: iKivelson et al.l Il997l : iThorne et al] 



I997I : iBespalov et all I2006I '). before being lost down the 



dusk flank of the magnetot ail v ia the pinching off of 
plasmoids (e.g. [Vasvliunaj 1 19831 : IWoch. Krupp fc Lagd 
l2002l : IVogt et al.ll2bl0l ) . 

As the plasma diffuses radially outward, its angular 
velocity drops (inversely with the square of the distance 
if no torques act) due to conservation of angular mo- 
mentum, such that a radial gradient of angular velocity 
is set up in the equatorial plasma. This angular velocity 




Figure 1. Sketch of a meridian cross section through a 
Jupiter-like exoplanet's inner and middle magnetosphere, show- 
ing the principal physical features involved. The arrowed solid 
lines indicate magnetic field lines, the arrowed dashed lines the 
magnetosphere-ionosphere coupling current system, and the dot- 
ted region the rotating disc of outflowing plasma. After Cowley 
& Bunce (2001). 



gradient, when mapped along the magnetic field to the 
ionosphere, causes an equatorward-directed (for Jupiter's 
magnetic field polarity) ionospheric Pedersen current to 
flow, the J X B force of which opposes the drag of the 
neutral atmosphere on the sub-rotating plasma. Angular 
momentum is transferred between the ionosphere and 
the equatorial plasma by the sweep-back of magnetic 
fleld lines into a lagging configuration, such that the 
ionospheric Pedersen current is balanced in the equatorial 
plane by an outward-directed (again, for Jupiter's magnetic 
field polarity) radial current, the J x B force of which 
tends to return the equatorial plasma back to corotation 
with the planet. Current continuity between these two 
field-perpendicular currents is maintained via field-aligned 
(Birkeland) currents, the inner upward component of which 
is thought to generate Jupiter's main auroral oval emission 
dCowlcv fc Buncel 12001) : Imiij I2OOII : ISouthwood fc KivelsonI 
[2001) . The current sys t em w as studied in detail theoretically 
by [Nichols fc CowlevI (|2003l ). who considered the effect of 
two poorly-constrained but important system parameters, 
the effective ionospheric Pedersen conductance Ep, and 
the plasma mass outflow rate M, and they derived analytic 
approximations appropriate for small and large radial 
distances, the former of which will be instrumental in the 
present work. INichols fc CowlevI went on to examine 

the effect on the current system of modulation of the 
ionospheric Peder sen conductance d u e to auror al electron 
precip itation, and INichols "fc" CowlevI (120051 ) and iRav et al.l 
(I2OI0I) have studied the e f fect o f field- aligned voltages. 
ICowlev. Nichob fc Andrewd (|2007l ) considered the effect on 
the current system of solar wind-induced expansions and 
contractions of Jupiter's magnetosphere, and recently the 
modulation of the current system by diurnal variation of 
the ionospheric Pedersen conduct ance caused by solar illu- 
mination has been considered by iTao. Fuiiwara fc Kasabal 
(2010) . 

In this paper we consider the application of the 
model describing Jupiter's magnetosphere-ionosphere cou- 
pling current system to Jupiter-like exoplanets with in- 
ternal plasma sources such as active moons. The stabil- 
ity of satellites about exoplanets ha s been considered (see 
e.g. iDomingos. Winter fc Yokovamal (|2006l ) and references 
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therein), since this is an important issue for the habitabil- 
ity of moons orbiting close-in 'hot Jupiters'. However, this 
is much less of a problem for satellites about more distant- 
orbiting exoplanets, which will be shown to be important for 
the present study. For example, iBarnes fc O'Bricnl (|2002l ) 
could not place any mass limits on satellites orbiting plan- 
ets beyond ~0.6 AU. The vulcanism of lo is driven by 
tidal dissipation, which also acts to dampen the eccentric- 
ity of the moon's orbit, such that the small eccentricities of 
the orbits of the Galilean satellite s of ~0 . 001-0. 01 are only 
raised due to Laplace resonance (jWeisd [2004i V The outer 
limit for stable orbits of a prograde moon was given by 
iDomingos, Winter &: Yokovamal l|2006h to be 



R„ 



X 0.4895(1 -1.0305ep-0.2738e™.) ,(1) 



where a is the semi-major axis of the planet's orbit, Mp 
is the mass of the planet. A/* is the mass of the parent 
star. Bp is the eccentricity of the planet's orbit and Cm 
is the eccentricity of the moon's orbit. For values appro- 
priate to lo (i.e. a = 5.2 AU, Mp = 1.899 x lO^'' kg, 
M* = Mq = 1.988 X 10^'' kg, ep = 0.048, and e„ = 0.004), 
_Rinax — 346 Rj, such that volcanic satellites orbiting 
jovian-like planets are expected to be stable (as is of course 
observed in the solar system, e.g. lo orbits at 5.9 Rj, 
well inside -Rmax). Indeed, for jovian mass planets orbiting 
beyond 1 AU, satellites with roughly lo's orbital distance 
will be typically stable, and only the planets with very 
eccentric orbits and high mass parent stars (e.g. a — 1 AU, 



20M, 



0' 



0.75, such that _R„ 



5.8 Rj) would 



be expected to be devoid of stable satellites at lo's orbital 
distance. 

Both giant planets in our solar system which have 
been extensively studied by in-situ spacecraft, i.e. Jupiter 
and Saturn, possess moons which actively outgas into the 
near-planetary space. At Jupiter this is lo, which emits 
material at the rate of ~1000 kg as discussed above, 
whilst at Saturn this is Enceladus, whose cryo-volcanoes 
emit water group ions at estimated rates ranging from a 
few to a few hund red kg s"^ (e.g. iKhurana et al.l I2OO7I : 
iPontius fc Hiill2006l ). It is thus reasonable to assume that 
such active moons may be relatively prevalent amongst 
Jupiter-like exoplanets, and the resulting implications for 
detectability of auroral radio emissions is therefore consid- 
ered here. We begin from the assumption that Jupiter-like 
exoplanets are strongly illuminated by their parent star, 
such that the ionospheric Pedersen conductance is high 
enough to maintain near-rigid corotation throughout the 
magnetosphere, a condition which, as reviewed below, 
maximises the field-aligned current density for a given mag- 
netosphere. We then compute the currents, and thus the 
resulting radio power output, for varying configurations of 
Jupiter-like exoplanets. Note that such strong illumination 
could be associated with close 'hot Jupiters', but equally 
applies to more distant planets orbiting active stars. The 
parameters examined are the plasma mass outflow rate, 
planetary orbital distance, and rotation rate. We show that, 
for planets with host stars more active than the Sun, only 
relatively modest modifications from the jovian system 
parameters are required to produce potentially-detectable 
configurations, and by doing so we open up the catalogue of 
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Figure 2. Plot showing (a) the magnitude of the north-south 
magnetic field strength threading the equatorial plane |-Bze| in 
nT, (b) the associated flux function Fe in nT Rp, and (c) the 
ionospheric co-latitude 6i to which the fleld lines threading the 
equatorial plane at radial distance pe map. The red solid lines 
show the 'baseline' model values given by Eqs. |4] and \5\ while 
the black long-dashed lines show the planetary dipole values for 
comparison. The other solid lines show the values modifled by us- 
ing Eqs. \E\ and |9] for assumed magnetopause standoff distances 
(Rmp/^p) ^ 85, and by simply extending Eqs. |4] and \5\ for 
(^imp/Rp) > 85. The colours blue, green, red, and black corre- 
spond to (i?mp/Rp) values of 20, 50, 85, and ^ 150, respectively. 
The horizontal dotted lines in panels (b) and (e) show the value 
of Foo ~ 2.841 X IC* nT Rp^, and the corresponding limiting 
ionospheric co-latitude of the outer boundary of the model at 
~ 15°. 



potential candidates for detection by radio telescopes such 
as LOFAR to a class of planet previously overlooked. 



2 THEORETICAL BACKGROUND 

In this section we discuss the equations which describe the 
system, along with the approximations appropriate to the 
case of Jupiter-like exoplanets. The a naly s is fol lo ws that 
which has been described previously bv lHilll (Il979l). Pontius 



(19971), ICowlev fc Buncd l|200ll ). and iNichols fc Cowlev 
(,2003i ), such that here we provide only a brief outline. 

2.1 'Baseline' magnetic field model 

We start by outlining the magnetic field model, which acts 
as the basis for the current system. It will become neces- 
sary to consider the variation of the current system with 
sub-solar magnetopause standoff distance Rmp, and here we 
thus describe the 'baseline' model which is then modified by 
changing Rmp- The magnetic field model is that which has 
been used previously for Jupiter bv lNichols fc CowlevI (|2004l . 



© 2011 RAS, MNRAS 000, [T]-?? 



4 J. D. Nichols 



l2005l ). It is an axisymmetric poloidal field with the magnetic 
axis co-aligned with the spin axis, such that it can be de- 
scribed by the flux function F{p, z), related to the magnetic 
field B by 



B 



(2) 



where p is the perpendicular distance from the magnetic 
axis, z is the distance along the axis from the magnetic equa- 
tor, and if is the azimuthal angle. Magnetic field lines are 
thus given by contours of F, such that mapping between the 
ionosphere (subscript i) and equator (subscript e) is simply 
achieved by writing Fi — Fe. For a near-dipolar ionospheric 
field we have 

Fi = -Be,/5? = BeqRl s\v? 9i , (3) 

where pi is the perpendicular distance from the magnetic 
axis in the ionosphere, B^q is the planet's equatorial dipole 
magnetic field strength (taken initially to be equal to 
Jupiter's value Bj — 426,40 nT in conformi t y wit h the 
VIP4 planetary field model of IConnernev et al.l (| 19981 )). Rp 
is the planetary radius (taken to be equal to Jupiter's radius 
Rj), and 6i is the magnetic co-latitude. In this model, the 
magnitude of the north-south magnetic field threading the 
equatorial plane is given by 



B 



exp 



,5/2 



+ B. 



°(^)' 



(4) 



where B'„ 



3.335 X 10^ nT, peo 



14.501 Ro 



Bo — 5.4 X 10"' nT, and m = 2. 71. This form clos e ly ap prox- 
imates th e field model used by Cowley & Bunco (2 00ll) and 
ICowlev. N ichols fc Bunco (2002); Cowley, Bunco fc Mchold 
l|2003h . who emplo yed the ' Voyager- 1 /Pioneer - 1 0' ('C AN') 
field model of IConnernev. Acufia fc Nes3 l|l98ll ) in 
the inner region and th e Voyager- 1 ('KK') model of 
[xii urana fc KivelsonI (|l993l ) in the outer region. This 'base- 
line' field model is shown by the solid red line in Fig. 
The values of \Bze\ are always less than those of the pure 
planetary dipole, shown by the dashed line in Fig. [2}3, due 
to the radial stretching of the field lines associated with 
the current sheet. The equatorial flux function is found by 
integration of Eq. [2] with the use of Eq. |4l such that 



f:{p.) = 

F, 



2.5pc 



(^) 



5/2 



(m-2) 



(5) 



where F^c ~ 2.841 x 10* nT Rp^ is the value at infinity, and 
r{a,z) = t°'~^e~* dt is the incomplete gamma function. 
This fiux function is shown by the solid red line in Fig. 
and is typically factors of ~5 greater than the corresponding 
dipole fiux function shown by the long-dashed line. 



2.2 Modification of the field structure by internal 
field strength and sub-solar magnetopause 
standoff distance 

We now consider how the above field model is modified 
by taking different values for the internal field strength 
and sub-solar magnetopause distance. As summarised by 



ICriefjmeier et all (|2004l ). there are a number of models which 
estimate the planetary field strength given the planetary 
rotation rate, the radius of the dynamo region, and the 
mass density an d conductivity within the dyn amo region, 
although recently iReiners fc Christensel] (|2010l ) used an ap- 
proach which was independent of the planetary rotation 
rate. Here we will assume for simplicity jovian values for 
the planetary radius, core conductivity and mass density, 
such that the planetary field strength is solely dependent 
on the planetary rotation rate Q.p. The planeta ry magnetic 
moment th us scale s as oc Op", where n = 1 jBussd 19761 : 
ISanolll993l). 3/4 (jMi^utani. Yamamoto fc Fuiimural Il992l) 



or 1/2 I Stevensor 19831 : iMizutani. Yamamoto fc Fuiimur; 
13). We thus take n = 3/4 as representative of these 



[199^. 

We thus take n 
choices, such that 

= (117) ' 



i 



and 



Fe = 



0.J 



f: 



(6) 



(7) 



where Q.j is Jupiter's planetary rotation period given by 
SI J = 1.76 X 10~* rad s"'^, and Bl^ and F* are given by 
Eqs. |4] and (5] respectively. 

Turning to the sub-solar magnetopause distance, this 
is expected to be dependent on the planet's orbital dis- 
tance from the star via the stellar wind dynamic pres- 



sure Ps 



as discussed below, and on the internal field 



strength. For the former, we e mploy the method of 
ICowlev. Nichols fc Andrewj (|2007l ). who studied the effect 
on the current system of changes to the magnetospheric size 
due to solar wind-induced compressions and expansions of 
Jupiter's magnetosphere. We note, however, that as well 
as considering the middle magnetosphere current sheet re- 
gion described by the above field model, those authors also 
incorporated 'outer magnetosphere' and open field line re- 
gions which are not included in this study since a signif- 
icant fiow shear at the open-closed field line boundary is 
not expected at pla nets w ith high ionospheri c con ductance, 
as considered here JlsbeU . Dossier fc Wait j Il984l ). Follow- 
ing ICwlexI^ichoir^&^^idrew^ l|2007l ). we set the reference 
boundary of the 'baseline' field model to be at 85 Rp, as 
shown by the red lines in Fig. [21 and modify the magnetic 
field model for smaller values of Fimp, by balancing the mag- 
netic flux lost due to the closer outer boundary with the 
addition of the equivalent flux via a uniform southward fleld 
AS throughout the system, such that 



■KRmp /^B = 2tv{F,{R„ 



F485 Rp)) 



(8) 



where AB thus represents the effect of the magnetopause 
currents as seen inside the magnetosphere. The fiux function 
F is then equivalently modified by the addition of a term 
AF given by 



(9) 



AF = F,(85 Rp) - F,(i?™p) + , 

such that an identical amount of fiux (Fe(85 Rp) ~ 
3.17 X 10'* nT Rp^) is contained within the boundary of the 
model for any value of Rmp smaller than 85 Rj, and thus 
the ionospheric location of the outer boundary is kept fixed 
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at ~ 15.8°. For values of Rmp larger than 85 Rp, we simply 
use Bze and Fe, since beyond ~130 Rp the magnitude of 
AB becomes larger than the more distant values of |-Bze|. 
Simply using Bze and Fe to larger radial distances does 
result in a modest increase in the amount of flux contained 
within the model, thus modifying the ionospheric location 
of the outer boundary, but it is obvious from the difference 
between the lowest values of the red lines and the horizontal 
dotted lines in Figs.[2}3 and c that this increase is very small, 
and the change in location of the open-closed field line 
boundary is less than 1°, i.e. less than obser ved changes in 
the l o cation of Jup i ter's main au roral oval l|Grodent et al.l 
I2OO8I : iNichols et atl 12001 l2009bl '). The magnetic field and 
flux values for these modified field structures are shown for 
magnetopause standoff distances of 20, 50, 85, and ^ 150 Rp 
by the blue, green, red, and black lines in Fig.[2l respectively. 

We now discuss how the magnetopause standoff dis- 
tance Rmp is expected to depend on planetary orbital dis- 
tance Rorb- The position of the magnetopause is dependent 
on the stellar wind dynamic pressure psw ~ Qsmvlw, where 
Vsw and Qsw are the stellar wind velocity and density, re- 
spectively. These parameters are related to the mass loss 
rate of the star, which is measurable by detecting Ly-a ab- 
sorption from the collision of the stellar wind w ith the sur- 
rounding interstellar medium (|Wood et al.ll2005l ) . These au- 
thors showed that stars more active than the Sun (which 
may be most important for internally-generated exoplane- 
tary radio emission, as discussed below) exhibit mass loss 
rates ranging between ~0. 01-100 times the solar value, and 
conjectured that although younger, more active stars are 
expected to exhibit greater mass loss rates, the high mag- 
netic field strengths of the most active stars may in fact 
act to inhibit stellar wind outflow. In view of these re- 
sults, we employ solar wind velocity and density values 
as reasonably representative of the large observed range 
of stellar mass loss rates for active stars. We thus assume 
that the stellar wind velocity Vsw is constant beyond 10 
stellar radii at 450 km s ~^, a value typica l of those ob- 
served at both ~1 AU (iHundhausenI Il995l ') and ~5 AU 
l|Nichols, Cowley fc McComasll2006l '). We then assume that 
the stellar wind density Qsw falls inversely with the square of 
the distance from the star (thus neglecting plasma sources 
in the interplanetary medium such as interstellar pick-up 
ions which become increasingly significant toward the outer 
regions of the heliosphere (e .g. Intriligato r. Siscoe fc Milleij 
[1996; McComas. Elliott fc S chwadron 201^)). The observed 
values of the solar wind density actually vary significantly 
over a solar rotation period due to the formation of coro- 
tating interaction regions in the interplanetary medium, 
and at Jupiter's orbit t akes values between 0.01 — 1 cm~"^ 
l|Nichols. Cowlev fc McComas 2006! ) . with typical values be- 
ing of order 0.1 cm~^. We thus employ the representative 
value of 0.1 cm~^ at 5.2 AU, such that the stellar wind dy- 
namic pressure varies with orbital distance as 



Psw 



0.9lfi-^VnPa 

V Rorb J 



(10) 



In the regime where the interplanetary magnetic field (IMF) 
strength is negligible with respect to the planetary field, 
stress balance between the dynamic pressure of the shocked 
stellar wind plasma in the magnetosheath adjacent to the 
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Figure 3. Plot showing (left axis) sub-solar magnetopause dis- 
tance Rmp and (right axis) stellar wind dynamic pressure versus 
orbital radius Roj-b computed using Eq. 1101 and the Huddleston 
et al. (1998) empirical relation for Jupiter given by Eg. 1131 (solid 
line) and the vacuum dipole relation given by Eq. [12] (dashed line), 
both for jovian values of the planetary magnetic moment. Also 
shown by the vertical dotted line at 5.2 AU is the range of magne- 
topause distances observed for Jupiter, thus the range over which 
the Huddleston et al. (1998) relation was defined, i.e. 45-100 Rj, 
and the horizontal dotted lines show the orbital radii to which 
this range corresponds in the present study, i.e. ~1.7-10 AU. 

magnetopause and the magnetic pressure of the compressed 
planetary vacuum dipole yields 



R 



Rp 
where 



swPsw 



1/6 



(11) 



» 0.88 for a high-Mach stellar wind 
( Spreiter fc Alksnel ll970') and km represents the factor by 
which the magnetopause field is enhanced by magnetopause 
currents, given by ~ 2.44 for a sub-solar boundary of re- 
alistic shape, i.e. lying between the values of 2 and 3 ap- 
prop riate to planar and spherical boun daries, respectively 
(e.g. iMead fc BeardI Il964l : Ulexeevl 120051 ') . Noting, however, 
that the precise values of ksm and km are in practise unim- 
portant due to the 1 /6 exponent in Eq. 1111 this relation is 
approximated by 



R-n 



Rp 



/^o psw 



1/6 



(12) 



similar relations to which have been used previously by 
variou s authors discussing ex o plane ta ry magneto s pheres 
(e.g. ICuntz. Saar fc MusielakI (20001: IZarka et all I2001I: 



Griefimeier et al.ll2004 l2005l: [Griefimeier. Zarka fc Spreeuwl 



20071 : Ijardine fc CameronI 120081 ). However, it is well known 



that this formula underestimates the size of Jupiter's mag- 
netosphere by a factor of ~2, since Jupiter's magnetopause 
is inflated to ~45-100 Rj, with a mean of ~75 Rj, by the 
presence of internal plasma through the co mbined actions 
of thermal press ure and centrifugal forc e l|Khurana et al.l 
2004). Therefore. iHuddleston et al.l (119981 ) presented an em- 
pirical relation for the variation of Jupiter's sub-solar mag- 
netopause standoff distance Rmp with Psw, given by 



' Rn 



35.5 



Rj J (ps™/nPa)0-2 



(13) 



Exoplanets with active moons would similarly be expected 
to contain signiflcant internal plasma, and here we thus em- 
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ploy a modification of Eg. 1131 taking into account tlie depen- 
dence on tlie internal field strength in Eq. 1121 to estimate 
the values of Rmv, i-e. 



\r7) 



35.5 



1/3 



(14) 



We show in Fig. [3] the variation of sub-solar magne- 
topause distance Rmp and stellar wind dynamic pressure 
Psw with orbital radius Ror t computed using Eq. [10] and 
the iHuddleston et al] (Il998l ) empirical relation for Jupiter 
given by Eq. [T3] (solid line), along with the vacuum dipole 
relation given by Eq. [12] (dashed line). The vertical dot- 
ted line at 5.2 AU shows the range of sub-solar magne- 
topause distances observed for Jupiter, i.e. ~45-100 Rj 
corresponding to a solar wind dynamic pressure range of 
~0. 01-0. 32 nPa. henc e represen ting the range over which 
the empirical .Huddlcston et al.l ((1998) relation is known to 
be valid. The horizontal dotted lines then show the orbital 
radii to which this dynamic pressure range corresponds in 
the present study through Eg. 1101 i.e. ~1.7-10 AU, although 
we acknowledge that intrinsic variations in the stellar wind 
dynamic pressure will actually induce a range in the sub- 
solar magnetopause standoff distances of an exoplanet at 
a given orbital distance, as observed at all solar system 
magnetospheres, an effect which we do not consider here. 
The validity of an empirical relation such as Eq. [T3] be- 
yond the observed values i s, of c ourse, not known. How- 
ever, the iHuddleston et all (|l998l ) relation encapsulates in 
a relatively simple manner the response of a rotationally- 
driven, plasma-filled magnetosphere to changes in solar 
wind dynamic pressure, and produces results consistent with 
the 3D MHD simulations of Ju piter's magnetosphere by 
lOgino. Walker fc KivelsonI (| 19981 ) and, more crucially, with 
in situ observations of the boundary location. On the other 
hand, the vacuum dipole relation given by Eg . 1 1 2 1 represents 
an approximation at odds with a key property of the systems 
studied here, i.e. the radial distension of the field through 
the actions of the centrifugal force and internal plasma pres- 
sure, and does not provide results consistent with in situ 
observations of the jovian magnetosphere for any solar wind 
dynamic pressure value. Thus, in the absence of any infor- 
mation about the size of extra-solar jovian-type magneto- 
spheres, we use the expression validated by observations of 
Jupiter as a reasonable estimation, even though this some- 
times involves extrapolation beyond the originally-observed 
range. 



2.3 Steady state plasma angular velocity and 
current system equations 

We next consider the angular velocity profiles arising from 
the application of Newton's second law to a radially- 
outward steady flow of plasma fro m th e torus source. 
As derived previously b y iHilll (1 19791) and jPontius 



and discussed further by ICowlev, Nichols fc Buncd (l2002f): 



1997 ). 



ICowlev. Nichols fc Andrewsl (|2007^ and [Nichols fc CowlevI 
(2003, 2004, 2005), the equatorial plasma angular velocity 
uj obeys the 'Hill-Pontius' differential equation 



Pe d 

2 dpe 



IB. 



M 



UJ 



(15) 



The left-hand side (LHS) of Eq. [15] represents the radial 
gradient of the plasma angular momentum per unit mass, 
whilst the right-hand side (RHS) is the ionospheric torque 
on the equatorial plasma per unit mass. We also note that 
the effective Pedersen conductance Ep is reduced from the 
true value Ep by Ep = (1 — fe)Ep, where the parameter 
k is related to the reduction of the angular velocity of the 
neutral at mosphere (ilt,) from rigid corotation (Q r>) via 
'slippage' ([Huang fc Hil]||l989l : [Millward et al.ll2005[ l. such 
that {flp — flp) — k{flp — uj). The value of k is somewhat 
uncertain, so in common with previous works we take 
k = 0.5. 

We now discuss the equations which describe the re- 
sulting magnetosphere-ionosphere coupling currents. First, 
taking the ionospheric field to be vertical and equal to 2Beq 
in strength (an approximation valid to within ~5% in our 
region of interest for planets with near-dipolar ionospheric 
fields (Nichols & Cowley 200>|)), the equatorward-directed 
height-integrated Pedersen current ip is given by 



(16) 



ip = 2Y,*pBeqO,ppi ( 1 — — 



Current continuity and the assumption of north-south sym- 
metry then yields for the equatorial radial current integrated 
across the width of the current sheet ip 

peip = 2piip . (17) 

From Egs. 1161 [17] and [3] and recalling that Fi = we have 



AT.*pFe9.p I , uo 



(18) 



such that the total radial current integrated in azimuth Ip 
is 



27rpeip 



: 87rE*pr2„Fe 1 



UJ 



(19) 



equal to twice the azimuth-integrated Pedersen current Ip = 
2n piip fiowing in each hemisphere. The field-aligned current 
density, e.g., northward out of the northern surface of the 
equatorial current sheet is then given, through current con- 
tinuity, by the divergence of the radial current, such that 



1 ^ . Id,., 
= -2^-*'' = -2^d^(^^^'') 



(20) 



where the factor of a half arises from the assumption of 
hemispheric symmetry, i.e. an equal and opposite current 
fiows out of the southern face of the current sheet. Noting 
that the quantity [j\\/B) is constant along a field line, such 
that {2Beq)) = (jze/Bze), we find the following expres- 
sion for the field-aligned current density at the top of the 
ionosphere 



Be 



dio 



27rpe|Bze| Ape 



(21) 



With these currents in mind, we now consider the approxi- 
mations appropriate to strongly-irradiated Jupiter-like exo- 
planets. 
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2.4 High-conductance approximation and 
conductance estimations 

Equation [15] is a first order, linear, inhomogeneous ordinary 
differential equation, which can be solved analyticall y for 
power law magnetic field models, such as a dipole field (|Hilll 
1 19791 . 1200 ih or the 'KK' field model valid throughout the ma- 
jority of Jupiter's middle magnetosphere (Nichols fc CowlevI 
I2OO3I ). and numerically for other field models. The solutions 
have the property that they are dependent on the quotient 
(Ep/M), as shown in Fig. [4^ for roughly jovian values of 
the ionospheric Pedersen conductance Ep and equatorial 
plasma mass outfiow rate M. Specifically, the typical scale 
over which the plasma falls from rigid corotation is called 
the 'Hill distance' ph after T. W. Hill, who first derived this 
scale length for a dipole magnetic field (not to be confused 
with the radius of the Hill sphere within which a body's 
gravitat ional field is dominant). The scale length was mod- 
ified bv lNichols fc CowlevI (|2003l ') for the case of a magnetic 
field that varies as a power law with arbitrary exponent m, 
and for which the field lines map from the equatorial plane 
to a narrow band in the ionosphere. This is appropriate for 
the stretched current sheet magnetic field of Jupiter's middle 
magnetosphere, as can be appreciated from the ionospheric 
mapping shown in Fig. [2J;, in which between 20-150 Rj the 
current sheet field models indicated by the solid lines map 
between ^15-17°, whereas the dipole field indicated by the 
dashed line maps to a much broader region between ~5-13°. 
The current sheet scale length is given by 



PH 



M 



l/m 



(22) 



where Fo is the value of the fiux function at the narrow iono- 
spheric band, taken here to be equal to Foo, shown by the 
dotted horizontal line in Fig. [S] Thus, for a given equatorial 
distance, the plasma an gular velocity incr e ases f or increas- 
ing values of (Ep/M). iNichols fc CowlevI (|2003l ) therefore 
derived analytic approximations to the solution appropriate 
to both the inner region, where the plasma near-rigidly coro- 
tates, and the outer region, where the plasma can be con- 
sidered to be either free of ionospheric torque or stagnant. 
The inner region appr oximation for the equator ial plasma 
angular velocity, which [Nichols fc CowlevI l|2003t ) suggested 
is valid out to a distance ps um of 



PS lii 



(23) 



is the leading term in the series solution of Eq. 1151 taking 
(M/Ep) as th e formal expansion param eter and (tj/ilp) = 1 
at pe = (see lNichols fc CowlevI (|2003l ) for further details), 
such that 



Pe_ 
pH 



(24) 



as shown by the dashed lines in Fig. Substitution of 
Eq. [24] into Eqs. [19] and [21] thus yields 



IpS ~ 2ip 5 — 



and 



J\\i S 



IB. 



TVpe\Bz 



dpe 



(25) 



(26) 



5x10"' 




Figure 4. Plot showing profiles of the current system parame- 
ters using the full solution of Eq. 1151 (solid lines) and the high- 
conductance approximation given by Eq. 1241 (dashed lines) . Pa- 
rameters shown are (a) the equatorial plasma angular velocity, 
(b) the azimuthally-integrated equatorial radial current, and (c) 
the field-aligned current density. Three results are shown for 
(SJ,/M) = 5 X 10-3, 10-3, and 2 x IQ-"' mho s kg-i, with M 
set to 1000 kg in panels b and c. Note only one dashed line 
is shown in panels b and c since the high-conductance approxi- 
mations (given by Eqs. l25l and l26l for the radial and field-aligned 
currents, respectively) are independent of Sp. The tick marks 
show the limit of validity of the approximation as suggested by 
Nichols % Cowley (2003). 



for the total radial current and ionospheric field-aligned 
current density, respectively. These limiting currents, 
which are dependent only on M and not Ep, are those 
required to maintain near-rigid corotation throughout the 
magnetosphere, and can thus be viewed alternatively as a 
high-Ep approximation for a given value of M. 

The currents Ip and jy^ derived from the angular veloc- 
ity profiles shown in Fig. [3^ are shown by the solid lines in 
Figs. [4]d and c, respectively, for roughly jovian field and cur- 
rent system parameters, i.e. M=1000 kg s"^ and Ep=0.2, 
1, and 5 mho. It is apparent that the current magnitudes 
increase with increasing Ep, asymptoting toward the lim- 
iting currents shown by the dashed lines. Thus, the high- 
conductance approximation represents the condition under 
which maximum currents, and therefore maximum output 
radio power, are generated. The minimum conductance re- 
quired such that the high-conductance approximation is 
valid throughout the magnetosphere is found by setting 



■ ps lim from Eq. 1231 such that 



Rm 
R„ 



0.516 



PH 

Rv 



(27) 
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which, upon substitution of Eg. I22l and rearranging for Ep, 
yields 



: 0.956 



Rrr 



M 



Bo Fo y Rp 



mho 



(28) 



Since we are interested in current systems powerful enough 
to be detected across interstellar distances, as possibly pos- 
sessed by Jupiter-like exoplanets strongly irradiated by their 
parent stars, such a high-conductance is clearly appropri- 
ate and will be used as the basis for the results which follow. 

Before considering the precipitating electron energy flux 
and associated radio power output, it is appropriate here 
to discuss the region of validity of the high-conductance 
approximation by estimating how the ionospheric Peder- 
sen conductance varies with planetary orbital distance Rorb, 
and the type of star. We use a model for the Pedersen con- 
ductance which comprises two components, i.e. conductance 
produced by solar illumination, and auroral particle precip- 
itation. For the former, in the absence of more complex 
ionospheric models for Jupiter-like exoplanets we assume 
that the Pedersen conductance is, like the Earth's, linearly 
proportional to the ionospheric electron number density, a 
quantity proportional to the square root of the electron pro- 
duction rate, which is in turn linearly proportional to the 
solar ext reme ultraviolet flux (iRobinson fc VondrakI [T984I : 
iLuhmarJ 1995), although we note that actually the Earth's 
upper atmosphere is ionised by solar radiation from X-ray 
to FUV wavelengths (a range denoted by the term 'XUV'). 
Thus, we have for the Pedersen conductance generated by 
stellar illumination 

1/2 



2.6 

Rorb 



mho 



(29) 



,Lxuv , 

where Lxuv * and L-^^jjy q are the XUV luminosities 
of the exoplanet's parent star and the Sun, respectively, 
and where the constant of proportionality is set such that 
XUV ~ 0.5 mho at 5.2 AU for solar X UV luminosity, 
consi stent with estimates for Jupiter (e.g. lAchilleos et al.l 
Il998l ). The Pedersen conductance is therefore expected to be 
higher for stars with higher XUV luminosities, such as more 
active stars. For example, the ~170 nm flux from F2V-type 
stars i s two orders of magnitude larger th an solar-type G2V 
stars (|Kasting. Whittet fc SheldonI Il996l ). since G2V type 
stars are sub-luminous in the FUV by a factor of t en or 
more due to absorption in this waveband (|Lean|[l997l '). Lu- 
minosity at X-ray-EUV wavelengths is mainly non-thermal 
coronal emission, and surveys indicate that most coronal 
sources in the EUV are active stars with coronal X-ray and 
EUV luminosities several orders of magnitu de greater than 
that of the Sun (|Barstow fc Holberd 12001 ) . SteUar age is 
an importa nt factor which d etermines stellar X-ray-EUV 
luminosity (|Ribas et al.ll2005l '). with younger main sequence 
stars typically exhibiting higher X-ray-EUV emissions. 
In the results which follow, we will thus consider stellar 
XUV luminosities from 1 to IQpLxnv a r ange consistent 
with values observed by iRibas et al.l (|2005h for solar-type 
stars, and thus could be conservative with respect to other 
spectral types. 

As well as stellar illumination, Pedersen conductance is 
expected to be generated by auroral electron precipitation. 
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Figure 5. Plot showing the field aligned current normalised 
to the limiting current produced by field-aligned voltage $ 
normalised to the unaccelerated population thermal energy W^^, 
for the three current- voltage relations given by Egs. 1331 for ijp = 
16, appropriate for an accelerator height of ~2.5 Rj (dashed line), 
[34] (solid line) , and 1351 (dotted line). "Rel." means "relativistic" , 
"Non-rel." means "non-relativistic" , and Rb is the ratio of the 
ionospheric magnetic field strength to the field strength at the 
top of the voltage drop. 



Computations bv lMillward et al] l|2002l ) indicated that up to 
^8 mho is generat e d by a uroral precipitation at Jupiter, and 
iNichols fc CowlevI (|2004l ) used these results to compute the 
expected variation of conductance with field-aligned current 
density assuming a representative primary and secondary 
precipitating electron spectrum. In their model the conduc- 
tance plateaus at ~3 mho at ~0.5 \xK m^^, since at higher 
currents the primary electrons precipitate lower than the 
Pedersen conducting layer and the conductance is generated 
by the lower energy secondaries. A full analysis of the effect 
of precipitation induced conductance requires that Eq. 1151 
a nd Eq. | 21 | are solved simultaneously, as in the computations 
of iNichols fc CowlevI { 20041 ) ■ but here we simply assert that 
is everywhere greater than 0.5 (J-A m~^, as is generally 
the case for the powerful exoplanetary systems we consider 
below, such that we assume a constant auroral contribution 
to the conductance of 3 mho, such that 



L 



XUV i 



^xuv 



1/2 



2.6 

Rorb 



+ 3 mho 



(30) 



The maximum orbital distance -R*^^ 



for which the high- 
conductance approximation is valid is thus obtained simply 
by setting Ep given by Eq. [28] equal to Ep given by 
Eq. [3S1 and solving for Rorb, numerically in practice. This 
then limits the maximum power output for a given system, 
since at larger distances the high-conductance approxima- 
tion is not valid and the magnetosphere will not rigidly 
corotate, and, as will be seen below, at smaller distances 
the power required to maintain the more compressed mag- 
netosphere rigidly rotating is reduced. 



2.5 Parallel voltage, precipitating electron energy 
flux, and radio power output 

We now briefly discuss the precipitating electron energy flux 
and associated radio power output that results from the 
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field-aligned current given by Eq. (26] The maximum field- 
aligned current that can be carried by unaccelerated precip- 
itating magnetospheric electrons is 



J||io 



eTV 



Wth 



1/2 



V 2'nme 

with a corresponding precipitating energy flux of 



/ Wth 

Efo = 2NWth TT-^ 



1/2 



(31) 



(32) 



where e is the elementary charge, A'^ is the magnetospheric 
electron number density, Wth is the electron thermal energy, 
and me is the electron mass. In the absence of detailed mod- 
els or empirical values for exoplanetary magnetospheric plas- 
mas we take, in common with previous works, jovian values 
for the plasma parameters outside the current sheet, i.e. A'^ = 
0.01 cm"^ and Wth = 2.5 keV (|Scudder. Sittler fc Bridgl 
ligSlI ). such that ^ 0.013 iJ.A m"^ A significant field- 
aligned voltage is thus required to drive field-aligned cur- 
rents greater than this value. Th e current-voltage relation 
was first derived bv lKnightl j 19731 ) and is given, for an arbi- 
trary value of the ratio Rb of the ionospheric magnetic field 
strength to the field strength at the top of the voltage drop, 
by 



^ ) = 1 + (i?B - 1) oxp 



WthiRB - 1) 



(33) 



where $ is the field-aligned voltage. In the limit that the 
field-aligned voltage region is located at high altitude (typ- 
ically > 4 Rj up the field lines at Jupiter), the ratio Rb 
becomes large (e.g. Rb = 216 at 6 Rj) such that Eq. 1331 
reduces to 



j\\io 



= 1 + 



the relation 
ll200ll) 



used 



Cowlev. Nichols fc Bunce 



ICowlev. Nichols fc Andrewj (|2007l ) 



(34) 

Cowlev fc Buncd 



II2002D. 

Nichols fc Cowlev! 



l|2003l . I2OO4T2OO5I '). On the other hand. iRav et all (|2010l ) 
suggested that the accelerator would be fixed at low 
altitude (2-3 Rj), such that th e y emp loyed Eq. [321 with 
Rb = 16. However, the iKnightl 1 19731 ) relations given by 
Eqs. a-nd [M] assume the accelerating voltage is small 
such that the precipitating particles are accelerated to 
non-relativistic speeds, an assumption which is not valid 
when the field-aligned voltage approaches or exceeds the 
electron rest energy (~511 keV), as would be the case 
for powerful exoplanetary systems, as s hown below. This 
restriction was removed bv lCowlevI (|2006l ). who showed that 
for a relativistic accelerated population the current-voltage 
relation becomes, in the high Rb limit. 



j\\io J \WthJ 



( W"f ) 



/ me \ 



+ 1 



(35) 



where c is the speed of light. The current- voltage relations 
given by Eqs. 1331 1341 and [35] are shown in Fig. [5] in 
which it is apparent that the relativistic and fixed-height 
accelerator cases produce competing deviations from the 
simple relation given by Eq. 1341 In light of the uncertainty 
surrounding the exact plasma parameters at Jupiter-like 
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Figure 6. Plot showing (on the left axis) the equatorial surface 
magnetic field strength Beq normalised to Jupiter's value Bj, and 
(on the right axis) the estimated radio emission bandwidth versus 
planetary angular velocity over the range 0.3 < Qp/Qj < 10, 
computed from Eqs. [Gland 1401 respectively. 



exoplanets, we thus employ the limiting iKnightl (|l973l ) 
relation given by Eq. [34] and shown by the solid line in 
Fig. [S] as providing a reasonable representation of the 
available choices for the current- voltage relation, and we 
simply note that this relation has been successfully used 
previously to compute auroral brightness and precipitating 
electron energy v alues for Jupiter that a re consistent with 
observations (e.g. ICowlev fc Bunce|[200ll ). 

The precipitating energy flux associated with the elec- 
trons acc elerated through voltage ^ as given by Eg. 1341 was 
shown bv lLundin fc Sandahj (j 19781 ) to be 



' 2 



-1-1 



2 



3\\ 



(36) 



where the approximation on the RHS is valid where 
^ J II io , a condition that is valid essentially everywhere 
in the models we investigate. The total precipitating power 
for each hemisphere is thus obtained by integration of Ef 
over the hemisphere, such that 



27ri?p sin^i Er dOi 



(37) 



Although the exact details of the mechanism by which ra- 
dio frequency waves are generated from unstable auroral 
electron distributions remain to be determined, the most 
likel y candidate is th e electron cyclotron maser instability 
fe.g. IWu fc Le3ll979l ). which has a generation efficiency of 
^1% (|Zarkalll998l ~ such that we take the total emitted radio 
power to be given by 



Pr = 



100 



(38) 



We finally assume that the radio emission is beamed into 
solid angles of 1.6 sr, in conformity with observations of the 
iovian HOM and DAM emissions (IZarka. Cecconi fc Kurthl 
2004*). The effect of beaming is to increase the flux at a 
given distance but decrease the chance of detection, and 
here we assume that the beam from only one hemisphere is 
observable at any one time. Assuming such beamed emission 
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Figure 7. Plot showing tiie current system parameters versus 
ionospiieric co-latitude, for magnetopause standoff distances of 
Rmp = 20, 50, 85, and 200 Rp (blue, green, red, and black lines 
respectively), and with flp = flj and M = Mj. The labels in 
panel (a) correspond to the values of R^rb which result in magne- 
tospheres of these sizes, according to Egs. llOl and lf 31 Parameters 
shown are (a) the azimuthally-integrated ionospheric Pedersen 
current Ip in A as given by Eq. 1251 (b) the field-aligned current 
density at the top of the ionosphere in A m~-^ given by Eg. 1261 
(c) the minimum field-aligned voltage 3>min in V required to drive 
the current in panel (b) as given by Eq. 1341 and (d) the precipi- 
tating electron energy flux Ef in W m~^ as given by Eg. 1361 The 
labels in panel (d) show the radio power obtained by integrating 
the energy fluxes using Eqs. 1371 and 1381 Note there are no blue 
lines in panels (c) and (d) since for this case the current density 
in panel (b) is everywhere below the minimum current jyio for 
which field-aligned voltages are required (shown by the horizontal 
dotted line in panel (b)). 



over the bandwidth Av, the spectral flux density at a given 
distance s from the source is 

Pr 



= 1 « 2 A ' (39) 
1.6s-' 

where the bandwidth Av is taken to be equal to the electron 
cyclotron frequency in the polar ionosphere, i.e. 

eBea 



Au = 



(40) 



again assuming a uniform polar ionospheric field strength 
equal to 2B^q. We show in Fig. [S] the field strengths and 
resulting estimated radio bandwidths for the range of plan- 
etary angular velocities considered in this paper (see below). 



3 RESULTS 

We now consider the application of the analysis presented in 
Sec.[2]to various configurations of exoplanets. We first show 
in Fig. [7] the effect on the current system of taking different 



values of Rmp, whilst maintaining jovian values for the 
planetary angular velocity Sip = Q.j (and thus Beq — Bj) 
and mass plasma outfiow rate M — Mj = 1000 kg s~^, 
thus simulating transplanting the jovian system to different 
orbital distances, assuming that the XUV flux is at all 
distances sufliciently high such that high-conductance 
approximation holds everywhere. It is apparent from Fig. [7] 
that for values of Rmp less than the reference value of 
85 Rj the current amplitudes are, as expected from e.g. 
Eq. 1251 reduced owing to the elevated magnetic fleld 
strength |i3ze|, which we recall is increased via Eq. |8]for a 
smaller magnetosphere due to fiux conservation. For these 
model parameters, very compressed magnetospheres with 
Rmp smaller than ~20 Rj, are such that the field-aligned 
current density is everywhere less than j\\io, hence no 
field-aligned voltages are required and there is no bright 
discrete auroral emission. Hence, even if we neglect tidal- 
locking which would in reality act to prevent the fast 
rotation of such planets (discussed further below), very 
close-orbiting (Rort ^ 0.25 AU) 'hot Jupiter' systems 
jovian-like in all respects apart from the orbital distance 
would exhibit no main auroral oval auroras. On the other 
hand, similar magnetospheres larger than 85 Rj exhibit 
common current profiles, extending to further distances and 
thus modestly further toward the pole, depending on the 
size of the magnetosphere. The monotonically-decreasing 
field strength with increasing equatorial radius results in an 
associated increasing of the current magnitudes toward the 
pole, and thus precipitating electron energy fiux, toward 
the outer boundary of the model. Increasing the size of the 
magnetosphere thus increases the emitted radio power, as 
shown in Fig. [7jl. 

Having briefiy examined the effect of magnetospheric 
size, we now extend the the analysis by considering the 
parameters M and flp. There are presently no published 
theoretical limitations on the maximum plasma mass out- 
fiow rate from exoplanetary moons, which presumably would 
vary with the number of volcanic moons in a given system as 
well as the orbital and geological parameters of each satel- 
lite and planet. The theory considered in this paper is valid 
for all non-zero values of the plasma mass outfiow rate, and 
it is worth noting that at Jupiter the ionosphere and so- 
lar wind also act as plasma sources with rates of < 10 ^ 
and > 20 kg s"\ respectively l|Hm. Dessler fc Goert3ll983l '). 
such that in practice the mass outfiow rate will probably be 
non-zero in most systems, even those devoid of active moons. 
It therefore seems reasonable to compute the radio power 
emitted for a range of M values within an order of mag- 
nitude of the canonical jovian figure of 1000 kgs~^. The 
rotation rate of a planet is governed by its initial angular 
velocity and subsequent tidal dissipation. The origin of the 
former is poorly understood, but candidates include the rel- 
ative motion of planetesimals during protoplanetary accre- 
tion and, for the gas giants specifically, h ydrodynamic fi ows 
during accretion of hydrogen and helium (|Lissaueij|l993l l. A 
crude estimation of the maximum angular velocity allowed 
by centrifugal stability of a planet is 



( GMp 



1/2 



(41) 
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(B) M = 3 Mj , Op = Q, 




(C) M = 0.5 , Op = 5 Q, 
o loot ^ ^ ' ~ 




(J) Af = 5 , Op = 5 




Figure 8. Plot showing results using four combinations of the 
input parameters Qp and M,i.e. (a) {{M/Mj) = 0.3, (Qp/Qj) = 
1), (b) m/Mj) = 3, (np/nj) = 1), (c) m/Mj) = 0.3, 
(Hp/Hj) = 3), and (d) {(M/Mj) = 3, (S^p/Hj) = 3). Each 
section consists of two panels. The top shows the limiting con- 
ductance given by Eq. 1281 (solid line) and the estimated Peder- 
sen conductance in mho given by Eq. 1301 for four values of the 
stellar XUV luminosity, i.e. {Lxuv * / Lxuv &) ~ IjlOilOO, and 
1000 (dashed lines from bottom to top, respectively). The bottom 
shows with the left axis the radio power in W emitted, obtained by 
integrating the precipitating electron energy fluxes using Eqs. 1371 
and 1381 and with the right axis the maximum distance in par- 
sees at which sources emitting such powers would be detectable, 
assuming a detection threshold of 1 mjy. All parameters are plot- 
ted versus orbital distance Rorb- The vertical dotted lines in each 
section indicate where the estimated Pedersen conductances equal 
the limiting conductance, i.e. show the maximum distances R* , 
to which the high-conductance approximation is valid, and thus 
the maximum radio powers for these parameters. The lines in the 
bottom of each panel are shown dotted beyond the maximum val- 
ues of R*^fj for {Lxuv -k/Lxuv &) ~ lOOOi indicating that this 
region corresponds to XUV luminosities above the upper limit for 
a young solar-type star as indicated by Ribas et al. (2005) . 



where G is the gravitational constant equal to ~ 6.67 x 
10~^^ N kg~^, and which for a planet of Jupiter's mass 
and radius yields ^Imax — 3.3 Qj. However, Jupiter lies 
roughly on the boundary between the regime of bodies sup- 
ported by Coulomb pressure, for which Mp oc Rp, and de- 
generacy pressure, for which Mp oc Rp^, such that Jupiter's 
radius is near maximal for bodies of solar composition. The 



net effect is that all bodies from solar composition giant 
planets through brown dwarfs to the very lowest mass stars 
are expected to have radii similar to Jupiter, and it is thus 
possible for planets with mass of, e.g. 10i\fj to have angu- 
lar velocities up to ~11 fij. The time scale Tayn required to 
de-spin a planet by tidal dissipation is given by 




where Q is the planet's tidal dissipation factor typically 
given by Q = 10^ for jovian planets, Q.orb is the Keplerian 
angular velocity of the planet, and w here we note that 
factors of orde r unity have been omitted ('Goldre ich fc Soterl 
1966, ; Showman fc GuillotI liooai ). For planets of jovian 
parameters orbiting solar-type stars at 0.1, 1, and 100 AU, 
this expression yields tidal synchronisation time scales of 
10*, 10^*, and 10^" y, respectively, while for planets more 
susceptible to tidal locking e.g. those with Alp — O.lMj and 
A/* — IOMq, the synchronisation time scales are modified 
to 10^, 10^^, and 10^^ y. Thus, close-orbiting 'hot Jupiters' 
are expected to rapidly become tidally locked, but for 
planets orbiting beyond 1 AU, tidal effects will be negligible 
over time scales of e.g. the present age of the solar system 
of ~ 4.6 X 10^ y, such that the planetary angular velocity 
may not be much reduced from its initial value at formation. 

With the above considerations in mind, we now examine 
the system using four spot combinations of the param- 
eters M and Qj, i.e. {{M/Mj) = 0.3, (fip/nj) = 1), 
((M/Mj) = 3, (np/Qj) = 1), {{M/Mj) = 0.3, 
{Qp/flj) = 3), and {(M/Mj) = 3, (^p/^j) = 3). 
Results are shown in Fig. |8] which is divided into four 
sections of two panels corresponding to the above pairs of 
parameters. The solid line in the top panel of each section 
shows the limiting conductance given by Eq. 1281 The 
dependence of the limiting conductance on the parameters 
Rorb, M, and Q.p, obtained from Eqs. H [3 [TSl [HI andESl 
is approximately Sp oc MQ,p~°'^^ Rorb^'^ , such that it 
increases slightly faster than linearly with distance in each 
panel, and compared to the values shown Fig. |5^, those 
in Figs. [HJd, c, and d are thus multiplied by factors of 10, 
~0.41, and '-^4.1, respectively. The estimated Pedersen con- 
ductance is also shown by the dashed lines in each panel for 
stellar XUV luminosities of {Lxuv -*/ Lxuv q) ~ 1,10,100, 
and 1000 (from bottom to top, respectively). Despite these 
XUV luminosities being orders of magnitude larger than 
the solar value, since Ep oc {Lxuv */ Lxuv q)^^"^ 
conductance values are multiplied from the solar values by 
factors of only 1, 3.2, 10, and 32, respectively. However, 
it is obvious that the increased conductance generated by 
the higher XUV luminosities greatly increases the maxi- 
mum orbital distance R^^-b which the high-conductance 
approximation holds, as shown in each section by the four 
vertical dotted lines, for example in section (a) at ~4.4, 5.3, 
7.3, and 10.9 AU corresponding to conductances of ^^3.6, 
4.6, 6.6, and 10.6 mho, respectively. As shown in Figs. [8]d, 
c, and d, these values of R^rb are then modified from, e.g., 
the high XUV luminosity case of ~10.9 AU to ~3.5, 17.7, 
and 5.4 AU, due to increased M, fip, and both M and 57p, 
respectively. 

The effects of these changes on the maximum radio 
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power are shown in the bottom panels of each section of 
Fig. [8] The left hand axes indicate the emitted powers, 
obtained by integrating the precipitating electron energy 
fluxes using Eas. l37l and l38l whilst the right hand axes show 
the distance s in parsecs obtained from Egs. 1391 and 1401 with 
these powers and an assumed spectral flux density of 1 mjy. 
Note that the right hand axes in Figs.[8j: and d are different 
to those in Figs. [8^ and b, since the spectral flux density 
is dependent on Av, which is proportional to flp^* through 
Bsq. As also shown in Fig. [T] the radio power increases 
with orbital distance, with a slight discontinuity in the 
gradient where Rmp = 85 Rp. This discontinuity is simply 
an artefact introduced by the change in the behaviour of the 
equatorial magnetic field model at the reference boundary 
of the 'baseline' magnetic field model. For Rmp < 85 Rp, as 
the radius of the magnetosphere increases with increasing 
orbital distance Rorb, the equatorial magnetic field strength 
decreases due to flux conservation, such that the current in- 
tensities and power output increase as discussed above (c.f. 
the blue and green lines in Fig. [7|. On the other hand, as 
the radius of the magnetosphere further increases such that 
Rmp > 85 Rp the current intensity at a given co-latitude 
remains the same, but extends increasingly toward the pole 
(c.f. the red and black lines in Fig. [7|. The latter has less 
of an effect that the former, such that the power output 
increases less quickly with Rorb after Rmp passes 85 Rp. 
Thus, for the values of M and ilp in Fig. [HK, the maximum 
radio powers P,- for the four XUV luminosity cases are 
~ 7.2 X 10^°, f.8 X 10", 8.2 x lO", and 2.5 x 10^^ W, 
corresponding to maximum observable distances s of ~0.4, 
0.7, 1.5, and 2.6 pc, respectively. However, again taking 
the high XUV luminosity case as an example, this power 
is modified as shown in Figs. [8]d, c, and d to --^ 2.6 x 10^^, 
8.6 X 10^^, and 3.4 x 10^^ W, corresponding to maximum 
observable distances of ~2.7, 10, and 20 pc, respectively. 
Thus, from these four spot combinations of A/ and Qp it 
is apparent that the maximum radio power Pr is relatively 
insensitive to changes in M, since the higher powers avail- 
able at a given orbital distance Rorb for increased M are 
compensated for by a decrease in the maximum distance to 
which the high-conductance approximation is valid Rlrb- 
On the other hand, whilst changing Qp does have competing 
effects on different components of the system, the overall 
behaviour is one in which increased flp results in increased 
power. We further note that, although, as discussed above, 
the compressed vacuum dipole relation for the sub-solar 
magnetopause standoff distance is not strictly appropriate 
for the rotationally-dri ven, plasma- filled magne tospheres 
we consider here, the iHuddleston et al] (|l998h relation 
was derived using observations of jovian values over the 
range ^45-100 Rj, and extrapolation beyond this range 
may lead to uncertainties. We have thus also computed 
the above values using Eq. [12] in place of Eq. 1141 and 
find that the results are qualitatively similar, except that 
the values of Rorbt Pr, and s are modified by factors of 
^2-4, ^0.1-0.4, and ^0. 3- 0.6, r especti vely. As noted by 
ICowlev. Nichols fc Bunc3 (|2002l ) and iNichols fc CowlevI 
l|2003h . the effect of the stretching of the planetary field 
from a dipolar configuration into a magnetodisc structure 
is to amplify the upward field-aligned current density 
associated with the aurora and radio emissions by 1-2 
orders of magnitude, such that the precipitating electron 
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Figure 9. Coloured plots indicating the maximum radio pow- 
ers Pr available using different pairs of system parameters M 
and f2p for, in panels {a)-(d) respectively, tlie XUV luminos- 
ity cases (Lxuv */ Lxuv q) ~ 1,10,100, and 1000. Also shown 
by the black contours are the maximum distances s in pc at 
which sources of these powers are observable, assuming a detec- 
tion threshold of 1 mJy, and the white contours show the orbital 
distances R* , in AU at which these maximum powers are avail- 

orb ^ 

able. Labelled crosses refer to the spot combinations shown in 
Fig.E] 



energy flux, a quantity proportional to the square of the 
field-aligned current density through Eq. 1361 is modified 
by 2-4 orders of magnitude. Using Eq. [12] in place of 
Eq. [14] compresses the equatorial magnetic field to a more 
dipole-like form, such that the radio power available at a 
given orbital radius is significantly reduced. However, the 
effect of a closer boundary is also to reduce the conductance 
required to maintain near-rigid corotation throughout the 
magnetosphere, such that the maximum power reduction is 
partly mitigated by increased orbital radius at which the 
maximum powers are available. However, we again note 
that Eq. [12] does not reproduce results consistent with 
observations of Jupiter's magnetosphere, such that in the 
dis cussion that follows we employ results obtained using 
the iHuddleston et all l| 19981 ') relation. 

We now extend the analysis to cover the full parameter 
space within the considerations discussed above. We thus ex- 
amine the behaviour of the system over 0.1 < (M/Mj) < 10 
and 0.3 < {Qp/Qj) < 10. We should note that angular ve- 
locities above ^3 flj would require planets more massive 
than Jupiter and may thus have larger dipole moments, al- 
though as discussed above for simplicity we do not take this 
into account here and will be the subject of future study. 
Results are shown in Fig.[9]for the four different XUV lumi- 
nosity cases. The colour indicates the maximum radio power 
Pr available for each pair of parameter values, as shown in 
Fig. [8] by the powers at the locations of the vertical dotted 
lines, and the black contours show the maximum observable 
distances of sources emitting such radio powers, again as- 
suming a spectral flux density threshold of 1 mJy. Note the 
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kinks are an artefact of the model indicating the change of 
behaviour where Rmp = 85 Rj as discussed above. The white 
contours indicate the orbital distances R^rb ^-t which these 
powers are available, which we note increase with decreasing 
M and, to a lesser extent with increasing Qp, as previously 
noted. This plot confirms the overall behaviour apparent in 
Fig. [S] i.e. that the radio power is essentially independent of 
M, increases with Q,j, and increases with stellar XUV lumi- 
nosity. Thus, for stars with solar XUV luminosity, planets 
with (rjp/fij) ~ 5 are required to produce radio emissions 
detectable from beyond ~10 pc, but this is reduced to e.g. 
{^Ip/Qj) ~ 2 for stars with {Lxuv */Lxuv q) ~ 1000 and 
(M/Mj) ~ 2. In all XUV luminosity cases, a significant 
number of parameter combinations within an order of mag- 
nitude of the jovian values are capable of producing emis- 
sions observable beyond 1 pc, in most cases requiring exo- 
planets orbiting at distances between ~1 and 50 AU. For the 
two higher XUV luminosity cases, parameter combinations 
within an order of magnitude of Jupiter's could generate 
emissions detectable beyond ~50 pc. 



4 SUMMARY AND DISCUSSION 

In this paper we have provided the first consideration of 
magnetosphere-ionosphere coupling at Jupiter-like exo- 
planets. We have estimated the radio power emitted by 
such systems under the condition of near-rigid corotation 
throughout (a condition which maximises the field-aligned 
currents and thus the radio power for a given magneto- 
sphere), in order to examine the behaviour of the best 
candidates for detecting internally-generated radio emission 
with next generation radio telescopes such as LOFAR. We 
have thus estimated for different stellar XUV luminosity 
cases the orbital distances within which the ionospheric 
Pedersen conductance would be high enough to maintain 
near-rigid corotation, and we have then considered the 
magnitudes of the large-scale magnetosphere-ionosphere 
currents flowing within the systems, and the resulting 
radio powers, at such distances. We have also examined 
the effects of two key system parameters, i.e. the planetary 
angular velocity and the plasma mass outflow rate. 

The key results of the study can be summarised as fol- 
lows: 

(i) The radio power emitted increases with increasing sys- 
tem size, and thus increases with orbital distance within the 
limit of validity of the high-conductance approximation. 

(ii) The limiting orbital distance, which defines the maxi- 
mum radio power available for a given set of system parame- 
ters, increases with stellar XUV luminosity and planetary ro- 
tation rate, and decreases with magnetospheric plasma mass 
outflow rate. 

(iii) The overall effect is that the radio power emitted 
increases with planetary rotation rate, but is essentially in- 
dependent of plasma mass outflow rate since the higher pow- 
ers available at a given orbital distance for increased plasma 
mass outflow rate are compensated for by a decrease in the 
maximum orbital distance to which the high-conductance 
approximation is valid. 

(iv) In all XUV luminosity cases studied, a significant 



number of parameter combinations within an order of mag- 
nitude of the jovian values are capable of producing emis- 
sions observable beyond 1 pc, in most cases requiring exo- 
planets orbiting at distances between ~1 and 50 AU. For the 
higher XUV luminosity cases the observable distances for jo- 
vian mass planets can reach ~20 pc, and massive, rapidly 
rotating planets could be detectable beyond ~50 pc. 

However, we should note here the limitations of the 
simple model used in this study. First, we have not consid- 
ered how the structure of the magnetic field changes with 
planetary angular velocity and plasma mass outflow rate. As 
noted above, Jupiter's magnetosphere is^ partly inflated by 
the centrifugal force of iogenic plasma (|CaudaI 1986,'). such 
that systems with higher angular velocity or plasma loading 
would be expected to be further inflated, although this 
would be somewhat mitigated by the associated increase 
of the planetary field strength in the case of the angular 
velocity. However, as discussed in Section O the effect of the 
stretching of the planetary field from a dipolar configuration 
into a magnetodisc structure is to amplify the upward field- 
aligned current density associated with the aurora and radio 
emissions by 1-2 orders of magnitude, and therefore the 
power values derived here may be viewed as lower limits in 
this regard. Second, we note that the field-aligned voltages 
obtained by the current-voltage relation used in this study 
may underestimate the true voltages if the location of the 
accelerator is fixed at a few planetary radii up the the field 
lines, such that these power values might again be under- 
estimates. In a related point, we note that simply defining 
the plasma angular velocity as we have done ignores the 
effects of the significant field-aligned voltages, which act to 
de-couple the equatorial and ionos p heric plasma an gular 
velocities (|Nichols fc CowlevI l2005l : iRav et al.l I2OI0I '). an 
issue which should be studied in further works. Further, 
we note that we have not considered at all the stellar 
wind interaction mediated by magnetic reconnection at the 
dayside magnetopause, which at Jupiter could be associated 



polar auroras observed ( 


Waite et al.ll200ll:|Pallier & PranEe 


2001 


; Grodent et al.l 


2003bl: Bunce. Cowlev & Yeoman 



20041; iNichols et al.ll2009al ll3). and would thus sporadically 
increase the power output. Nor have we considered the 
effect of changing stellar wind dynamic pressure, which 
at Jupiter is known to modulate the in tensity of the 
UV and radio emissions by factors of ~3 (iGurnett et al] 
20021: ICowIev. Nichols fc Andrews! l2007l : IClarke et al.ll2009l : 
Nichols et al .Il2009bl '). and again, although the details of how 
this affects Jupiter's auroral and radio emissions remain to 
be fully determined, this effect may also act to increase the 
output powers from those computed here. However, we also 
reiterate that our assumption of solar values for the stellar 
wind velocity and density may represent underestimates for 
stars more active than the Sun, such as those much younger. 
Higher stellar wind dynamic pressure values would decrease 
the size of the magnetosphere for a given orbital distance, 
thus reducing the emitted radio power. It is clear that the 
evolution of internally-generated radio emissions over the 
lifetime of a star should therefore be considered in future 
studies. We should also recall that we have only considered 
here planetary magnetic fields with the same polarity as 
those of Jupiter and Saturn, i.e. with magnetic and spin 
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axes co-aligned to first order, since this is the configuration 
that has been most studied for bodies in our solar system. 
However, it is probable that only 50% of Jupiter-like 
exoplanets exhibit this polarity. The effect of reversed 
polarity is to reverse the direction of the current system, 
such that downward currents are replaced with upward 
currents and vice versa, such that the source population for 
the field-aligned currents in this case may be very different 
to tha t considered here (see e.g. iBunce. Cowley fc Yeomiml 
l|2004l for a discussion of field-aligned currents induced 
at Jupiter's dayside magnetopause) . The emissions from 
planets with the opposite polarity should be studied in 
future works. 

Finally, we note that the implication of these results 
is that the best candidates for detection of such internally- 
generated radio emissions are rapidly rotating Jupiter-like 
exoplanets orbiting stars with high XUV luminosity at or- 
bital distances beyond ~1 AU. This type of exoplanet has 
not previously been considered as potential detection can- 
didates for next generation radio telescopes, but searching 
for such emissions such may offer a new method of detec- 
tion of more distant-orbiting exoplanets less likely to be de- 
tected by those techniques which are biased toward close- 
orbiting 'hot Jupiters'. However, dual detections with radio 
telescopes and conventional methods would best constrain 
the planetary parameters. 
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